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Abstract 

We study FRW cosmology for scalar tensor theory where two scalar functions nonminimally 
coupled to the geometry and matter lagrangian. In a framework to study stability and attractor 
solutions of the model in the phase space, we simultaneously solve the dynamical system and 
best fit the stability and model parameters with the observational data. The approach imposes 
restrictions on the model constraints while providing information about the universe dynamics. 
The model predict current accelerating universe, with a phantom crossing in near future. 
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1. INTRODUCTION 



Recently, the observations of high redshift type la supernovae, the surveys of clusters 
of galaxies PP-[1], Sloan digital sky survey (SDSS) [S] and Chandra X-ray observatory [B] 
reveal the universe accelerating expansion. Also the observations of Cosmic Microwave 
Background (CMB) anisotropies [7j indicate that the universe is fiat and the total energy 
density is very close to the critical one p]. The observations though determines basic cos- 
mological parameters with high precisions and strongly indicates that the universe presently 
is dominated by a smoothly distributed and slowly varying dark energy (DE) component. 
A dynamical equation of state ( EoS) parameter that is connected directly to the evolution 
of the energy density in the universe and indirectly to the expansion of the Universe can be 
regarded as a suitable parameter to explain the universe acceleration [9]-[T2]. 

Motivated from string theories, the scalar-tensor models provide the simplest model- 
independent description of unification theories which predict couplings between scalars and 
curvature. They have assumed a prominent role since any unification scheme, such as 
supergravity, in the weak energy limit, or cosmological models of inflation such as chaotic 
inflation, seem to be supported by them [13]. In addition, they have been employed to study 
the current acceleration of the universe [Ti] |21j. 

On the other hand, by using the well-known geometric variables, Hubble parameter and 
deceleration parameter together with the new geometrical variables, the cosmological diag- 
nostic pair {r, s} ( or statefinder parameters) |22] . the acceleration expansion of the universe 
and differentiation among the cosmological models can be explained in order to better fit 
the observational data. The importance of the statefinder parameters to distinct DE cos- 
mological models is best realized, in particular, when considering the increased accuracy of 
the observational data during the last few years and generality of the DE models. These 
parameters, in a natural next step beyond the well known geometric variables, are to dif- 
ferentiate the expansion dynamics with higher derivatives of the scale factor and to explore 
a series of DE cosmological models, including A cold dark matter (LCDM), quintessence, 
coupled quintessence, Chaplygin gas, holographic dark energy models, braneworld models, 
and so on [23]-[26]. Moreover, since the cosmic acceleration affects the expansion history 
of the universe, to understand the true nature of the driving force, mapping of the cosmic 
expansion of the universe is very crucial [27j. Hence, one requires various observational 
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probes in different redshift ranges to understand the expansion history of the universe. One 
of these tests is the difference in distance modulus measurement of type la supernovae that 
helps us to testify the cosmological models. 

In this paper, we study the detailed dynamics of the nonminimally couples scalar field 
to gravity and matter lagrangian. Since the major difficulty in cosmological models is the 
nonlinearity of the field equations and thus limitation in obtaining the exact solutions, we 
investigate the asymptotic behaviour of the model, which provides the relevant features to 
be compared with the current physical data available for the universe. In this context, the 
perturbation methods, especially linear stability analysis which have been used to study the 
qualitative analysis of the equations and of the long term behavior of the solutions are being 
proposed in this work |28j |29j. 

The paper organized as follows: Sec. two is devoted to a detailed formulation of the 
model. In Sec. three, we obtain the autonomous equations of the model and the late time 
attractor solutions by using the phase plane analysis in addition to best fitting the model 
parameters. We test the model with observational analysis; examining the behavior of the 
effective EoS parameter of the model and statefinders. We also reconstruct the dynamics 
of the coupling functions to the curvature and matter lagrangian. In Sec. four, we present 
summary and remarks. 



2. NONMINIMALLY COUPLED SCALAR FIELD MODEL 

A general action in four dimensions where gravity and matter are nonminimally coupled 
to two different scalar functions given by 

A = j [G{<P)R - ]^<P,,<jy'^ + L^f{cp)]^gd'x, (1) 

where the four functions G'(0) and /(0) are not specified. Note that the four function G 
coupled to the curvature is different from gravitation coupling constant that appears in 
the standard general relativity. The action generalizes those used to construct extended 
inflationary models [30]. It is also conformally equivalent to the Brans-Dicke action in 
the Jordan frame by using the transformation if = G{(f)) and u){ip) = 2(dG{!fy]d<j))^ ^^"^ 
putting / = const. In our units, the standard Newton coupling is also recovered in the limit 
—\ [13j. The fleld equations can be derived by varying the action with respect to 
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where 



In addition, the variation with respect to gives the klein-Gordan equations 

In FRW cosmology, the Euler-lagrange equation, corresponding to the cosmological Einstein 
equations are and spatially homogeneous scalar field (f){t) are 



3i.^^-3^ + -, + --p (5) 

2H + 3H (6) 

where eq (|5| is the energy constraint corresponding to the (0,0)-Enistein equation. In addi- 
tion, by using the argument in [31j, for a spatially homogeneous scalar field 4>(t), the field 
equation is 

4> + 3H^ = 6{H + 2H^)^ - (1 - 37)p^/'(0). (7) 

where prime means derivative with respect to (f). In the above, we assumed that the universe 
is filled with the barotropic fiuid with the EoS to be Pm = 'ypm- From equations (5), (6) and 
(7), one can easily arrive at the modified conservation equation, 

(pj) + 3H{1 + -i)p,^f = (1 - 37)p„./> (8) 

In the following we assume that the matter presented in the universe is cold dark matter 
with 7 = 0. 



3. STABILITY ANALYSIS AND BEST FITTING 

In this section, we study the structure of the dynamical system via phase plane analysis, 
by introducing the following dimensionless variables, 

OH' 12GH^' 6GH^' 
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Using equations ([5])-([7]),and power law functions for G(0) = exp(a0) and /(0) = exp(/50) 
the dynamical equations in terms of the new variables become, 

r = r( A _ ^ _ A) (11) 

Z' = -3Z-XZ- 2Z{^) (12) 
where prime " ' " means derivative with respect to = ln{a) and we have, 

^ = (4y^)(-3 + -^- + -'f^-f^-31''+lJ^). (13) 

X , I3ZX , , 

^ xA+x^ (15) 



The Fridmann constraint equation ([s]) also becomes 



Z = l + X-F^ (16) 



Using the constraint (16), the equations (10)-(12) reduce to two coupled differential equa- 



tions for X and Y . Next, we obtain the critical points (fixed points) and study the stability 
of these points. Critical points are always exact constant solutions in the context of au- 
tonomous dynamical systems. These points are often the extreme points of the orbits and 
therefore describe the asymptotic behavior of the system. In the following, we find fixed 
points by simultaneously solving X' = and Y' = 0. Substituting linear perturbations 
X' X' + 6X', y — )■ y + 6Y' about the critical points into the two independent equations 
for X and Y, to the first orders in the perturbations, yield two eigenvalues Xi{i = 1,2). 
Stability requires the real part of all the eigenvalues to be negative. 

In the following we solve the reduced form of the coupled system of differential equations 
to find fixed points. From the numerical calculation we find that both critical points and 
eigenvalues in our model, depending on the model parameter a and /3, are highly nonlinear. 
In addition, the expressions for them are long and cumbersome such that we can not evaluate 
under what conditions the critical points are stable or unstable. In a different approach to 
the problem we solve the equations by simultaneously best fitting the stability parameters 



and initial conditions with the observational data using the method. This affords us to 
find solutions for the above equations and conditions for the stability of the critical points 
that are physically more meaningful and observationally more favored. In the next section, 
we solve the equations by simultaneously best fitting the model with the observational data 
for distance modulus. 

3.1. Observational best fitting with distance modulus, //(z) 

The difference between the absolute and apparent luminosity of a distance object is given 
by, fj,{z) = 25 + Slogj^o dL^z) where the Luminosity distance quantity, diiz) is given by 

In our model, from numerical computation one can obtain H{z) which can be used to 
evaluate fi{z). To best fit the model for the parameters a ,/3 the initial conditions X{0), 
Y{0), H{0) with the most recent observational data, the Type la supernovea (SNe la), we 
employe the method. We constrain the parameters including the initial conditions by 
minimizing the function given as 

^2 R Y(a\ v(^\ n(^\\ \- AX(0),y(0),g(0)) -/if-]2 

i=l « 



where the sum is over the SNe la sample. In relation (18), /i*'*^ and are the distance 
modulus parameters obtained from our model and from observation, respectively, and a is 
the estimated error of the ^f"^.Yiom. numerical computation. Table I shows the best best- 
fitted model parameters. 

TABLE I: Best-fitted model parameters and initial conditions, 
parameters a (5 X(0) y(0) i7(0) xl 



/2 . 



-0.23 -3.2 4.36 0.5 0.5 545.9688383 



The contour diagrams at the 68.3%, 95.4% and 99.7% confidence levels are given in 
FIG. 1). From the graph, one conclude that with 68.3%, 95.4% and 99.7% confidence level 
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the true values for both a and /3 he within the green, blue and red contours, respectively. 
Alternatively, we can plot the likelihood for the pair model parameters as shown in Fig. 2). 




Fig. 1: The constraint on a and /3 at the 68.3%, 95.4% and 99.7% confidence 

levels from Sne la for the model. 




Fig. 2: 1-dim and 2-dim likelihood for the model parameters 



In Fig. 3, the distance modulus, yu(-2), in our model is compared with the observational 
data for the obtained parameters and initial conditions using method. As can be seen the 
best fitted parameters and initial conditions are in good agreement with the observational 
data. 




0.5 1 1.5 

z 



Fig. 3: The best fitted distance modulus [i{z) plotted as function of redshift 
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In the following we investigate the stability of the model with respect to the best fitted 
model parameter. 

3.2. Stability of the critical points and phase space 

Solving the stability equations for the best fitted model parameter a and /3 we find eight 
fixed points with the stability properties as illustrated in tables II. As can be seen the critical 
points in the phase space are symmetric with respect to the symmetric axis X — This is 
due to introduction of the dimensionless variable which creates symmetric roots in the 
phase space. The fixed points FP7&8 lie on the symmetric line, thus, is a double fixed point. 

TABLE II: est fitted critical points 
points FPl FP2 FP3 FP4 FP5 FP6 FP7&:8 



X 0.5 0.5 0.9 0.9 1.5 1.5 

Y 0.6 -0.6 1.5 -1.5 0.9 -0.9 

property saddle point saddle point unstable unstable stable stable unstable 



Prom the above table we see that two of the critical points, FPl and FP2 are stable and 
the rest are unstable or saddle points. In Fig. 4) the trajectories leaving the unstable critical 
points FP5, FP6, FP7 and FP8, passing the saddle points FPl and FP2 and finally entering 
the stable critical points FP5 and FP6 in the phase plane are shown. 
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Fig. 4: The attractor property of the dynamical system in the phase plane. The trajectories 
with the model parameters a and /3 are shown by blue lines. 
The red trajectories correspond to the 
best fitted stability parameter and initial conditions. 

In Fig. 5, in more details for the given initial condition of F > 0, the trajectories (blue 
curves) approaching the stable critical point FP5 in the phase plane are shown. For the 
initial conditions of y < 0, then the trajectory shown entering the stable critical point 
FP6. The best fitted trajectories ( red curve) with the properties given in the previous 
section enter the stable critical points FP5 and FP6. These trajectories both fit the model 
parameters a and /3 and the initial conditions for X, Y and H. 
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Fig. 5: The attractor property of the dynamical system in the phase plane. The trajectories 
entering the best fitted stable critical points with the model parameters a and /3 
are shown by blue lines. The red trajectory approaching the FP5 and FP6 
corresponds to the best fitted stability parameter and initial conditions. 

COSMOLOGICAL PARAMETERS 



In order to understand the behavior of the universe and its dynamics we need to study the 
cosmological parameters such as EoS parameter. We have already verified our model with the 
current observational data via the distance modulus test. The EoS parameters analytically 
and/or numerically have been investigated by many authors for variety of cosmological 
models. Applying stability analysis and simultaneously best fitting the model with the 
observational data using method gives us a better understanding of the critical points. 
In our model, the effective EoS parameter is defined by Wg// — ~^~\^ where ^ is given 



in terms of the dynamical variables in equation (13). Among cosmological parameters, the 



statefinders are given by r = H / — 3g — 2 and s = (r — l)/3(q' — \) discussed here, where 
g in r and s is the deceleration parameter and in r in terms of new dynamical variables 
can be obtained by taking derivative of H. 

The calculated effective EoS parameter for the best fitted stability parameters a and /3 
and initial conditions exhibits phantom crossing behavior in the future where the universe 
is in an unstable saddle point FPl (Fig. 6). The current effective EoS parameter for the 
best fitted model parameters is cje// — —0.667. The graph also shows that in the far past 
for high redshift the universe is in matter dominated era in an unstable saddle point FP7 
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where Uefj = 0. The best fitted effective EoS parameter continues its journey towards 
unstable saddle point FPl at about redshift z ^ —0.3 where u^ff = — 1 in near future. The 
interesting point is that the unstable saddle point state of the universe is when it crosses 
the phantom divide line. Finally, the trajectory approaches the stable critical point FP5 in 
the future. 
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Fig. 6: The best fit graph of effective EoS parameter plotted as function of redshift 



Following in Chevallier-Polarski-Linder (CPL) parametrization model one can use 
linearly approximated EoS parameter, 



1 + ^' 



(19) 



where the parametrization is fitted for CPL parameters and Wq. In the following, using the 
method, we can best fit the CPL parameters with the observational data that have been 
used to find effective EoS parameter. Fig 7 shows the 68.3%, 95.4% and 99.7% confidence 
level for and uja that lie within the green, blue and red contours, respectively. 




0.14 0,34 0.54 0.74 0.94 1.14 1.34 1.54 



Fig. 7: The constraint on and uJa at the 68.3%, 95.4% and 99.7% confidence levels. 

Alternatively, the 1-dim and 2-dim likelihood for the CPL parameters and uja are 
plotted in Fig. 8. 
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Fig. 8: 1-dim and 2-dim likelihood for the CPL model parameters loq and uja- 

The best fitting procedure thus finds that the approximate values of the CPL parameters 
are uq = —0.72, Ua = 0.74. A comparison of the standard CPL parametrisation matching 
the best fitted effective EoS parameter is shown in Fig. 9. 
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Fig. 9: The best fitted effective EoS parameter in comparison with the CPL EoS 
parameterization with the best approximated parameters Uq = —0.72, Ua = 0.74 

By applying the statefinder diagnosis to the model. Fig 10 shows the best-fitted trajecto- 
ries of the statefinder diagrams {r, q}, {s, q} and {r, s}. From the graph it can be seen that 
the best-fitted trajectory passes LCDM state with {r,s} = {1,0} sometimes in the past. 
The current value of the best fitted trajectory and its location with respect to the LCDM 
state can also be observed in the {r, s} diagram. From Fig 10, we see that the current 
value of the statefinder {r, s} is close to the LCDM state and compatible with the recent 
observational data. 
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Fig. 10: The best fit graph of statefinders {r, g}, and {r,s} 

In Fig 11 we depict the corresponding dynamical behavior of the satefinder {r, s} against 
= ln{a). From Figs. 10 and 11 we observe that the universe starts its journey from 
unstable state in the past, passed the current state and eventually reaches a stable state 
in the future. Interestingly, the trajectories of the statefinders and their evolution show 
that sometimes in the future the universe approaches an unstable saddle point that is the 
extreme points of the statefinder trajectories corresponding to the state that the universe 
crosses phantom divide line. 
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Fig. 11: The best fit graph of statefinder parameters r and s plotted as function of = ln{a) 

We also reconstructed the best-fitted functions G{(j)) and f{<f)). From Fig 12, we see that 
the trajectories for the best fitted model parameters illustrate a monotonic increasing and 
decreasing behavior with respect to the scalar field (p for and G{(j)) respectively. 
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Fig. 12: The best fit reconstructed graph of G'(0)and plotted as function of 

5. SUMMARY AND REMARKS 

This paper is designed to study dynamics of the nonminimally coupled scalar field theory 
by using the stability analysis and the 2- dimensional phase space of the theory. In a different 
approach in stability analysis, here, we solve the system of autonomous differential equations 
by best fitting the model parameters and also the initial conditions with the observational 
data for distance modulus. Therefore all the critical points with the stability conditions 
are presented in the model are physically reliable and observationally verified. By stability 
analysis, we find two stable critical points in the model as shown in Fig. 4 and 5. With the 
best fitting results the trajectories approaching the stable points are observationally verified. 

We then study the cosmological parameters such as effective EoS parameter, u^jf for the 
model in terms of the dynamical variables introduced in the stability section. The result show 
that with best fitted stability parameter, the model exhibits an accelerating universe with a 
transition from We// > — 1/3 (decelerating universe) towg// < — 1/3 (accelerating universe)in 
the past. It also shows that the universe sometimes in the future passes an unstable state 
( saddle critical point in the phase space) where the effective EoS parameter crosses the 
phantom divide line. The best fitted effective EoS parameter is also compared with the 
CPL parametrisation of EoS and the best approximate value of the CPL parameters found 
to be uq = —0.72, Ua = 0.74. With the statefinder diagnostic, the best fitted statefinder 
parameters show that the current state of the universe is very close to LCDM. The statefinder 
graphs and also their evolution explicitly illustrate that the unstable saddle point is the 
extreme point of the statefinder trajectories. Finally, the best fitted reconstructed coupled 
scalar field functions, (p, G{(p) and /(0) against the scalar field display a smoothly decreasing 
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and increasing behavior. 
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